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984 01:21:28 VAX/VMS Macro V04-00 Page 
P- 38 on $}; 746 CMTHRTL.SRCIMTHDSINH.MAR; 1 . 
-TITLE MTHSOSINH 3 povene,! loating Hyperbolic Sine routine 
IDENT /1-008/ : File: MTHDSINH.MAR Edit: RNH1008 


s ebesendeseneesseserorbeonnnnsnesonsbensenaneorebensnnesenesensacenseseoseess 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


*® y 
‘@ ® 
*® ® 
*® ® 
‘® ® 
is THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
* ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE ® 
* INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
* COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
:® OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
:® TRANSFERRED. * 
* ® 
:® THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
;* AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT ®* 
:® CORPORATION. * 
:@ ® 
:® ® 
:® ® 
*® * 
:* * 
* * 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


LARA RARAAAAASASAAARAAAALEAARELASELAESEE RSE RA REESE EERE RARER SESE REESE RE RES RE RS SO 


; FACILITY: MATH LIBRARY 
; : "ABSTRACT: 


: S MTHSDSINH is a function which returns the double floating hyperbolic sine 
; of its double precision floating point argument. The call is standard 
; ; call-by-reference. 


; ; VERSION: 01 


: ; HISTORY: 
Peter Yuo, 29-Jun-77: Version 01 


: ; MODIFIED BY: 


“—— 
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-SBTTL HISTORY ; Detailed Current Edit History 


ALGORITHMIC DIFFERENCES FROM FP-11/C ROUTINE: none 
Edit History for Version 01 of MTHSDSINH 


0-2 he aa ere oe A * Silerecntpanes 
MTH$_... change 


v04-00 
2:46 (CMTHRTL. sR CIM THDSINH.MAR; 1 


Changed error endl Sea | mechanism. Put error result in RO:R1 before 


call 


Oo 
t 
is 


Six term Taylor series, in powers of argquaent , replaced 
uy eh x term Chebyshev series, in powers of ARG** 

th gverhang provided. This improves Ti om 18-May-1978; 
Mary 


Fix Loat3s constant. TNH 16-June-78 
Update version number and c rnk ee ht notice. JBS_ 16-NOvV-78 

- ¢ ange MTH__FLOOVEMAT to MTH Pioneer: JBS 07-DEC-78 

- Remove SSRADEF - not needed. “OSS -78 

- Add "*."' to the PSECT n't tives. nies -DEC-78 

- Declare externgis. 6, 

- Use MTHSDEXP_R6. SBL 7 Sep t- 15999 


IO 


ng MTHSS$SIGNAL in order to allow user modify error result. 


Changed lower Limit for Chebyshev approximation from 2**-27 to 


- Eliminated s second call to EXP for input values between 28.5*|n2 


and 1 


- Changed all final floating point divisions by 2 to interger 


subrtracts of 1 from the exponent field. 
- Extended maximum range from 87.69 to 128*ln2=88.72. 
- ite tne te ic for computing EXP(!x!=ln2) to reduce error. 


1-008 - Change W* to G* on call to MTH$$SIGNAL RNH 09-Sept-1981 


oe 
- 
ae 
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arative Part of Modul 6=-SEP 2:46 CMTHRTL.SRCIMTHDSINH.MAR; 1 
-SBTTL DECLARATIONS ; Declarative Part of Module 


; INCLUDE FILES: MTHJACKET.MAR 


EXTERNAL SYMBOLS: 


-DSABL GBL 

.EXTRN MTHSDEXP_R6 
EXTRN MTHSK_FLOOVEMAT 
~EXTRN MTHSSSIGNAL 


; EQUATED SYMBOLS: 


00004080 $D_1.0 *F1.0 3 1.0 
00000004 value 4 3 value.rd.r 
MACROS: none 


PSECT DECLARATIONS: 


-PSECT _MTHSCODE PIC,SHR,LONG,EXE,NOWRT 
; program section for math routines 


SOOoOQoQooooQoooooooooooooooooooo 


SOooooooooooooooo 
SOooOCCOoCoCoQooQooooooooooooooooo 
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So 

So 

So 
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; OWN STORAGE: none 


FP WWININIWINIWIIDINIPINININININIDNINIDN 2 9 2 POO OCOOODOOCOCOOOOOOOOOOOdMADCDODODOD -- 


1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
000 1 
000 1 
000 1 
000 1 
000 1 
000 1 CONSTANTS: 
000 1 
000 1 
000 1 D_127_L0G_2: 
2BDAC7E2 0F 334380 000 1 ~QUA *X2BDAC7E20F 3343B0 3 127*Ln2 
B28 121 D_128_L0G_2: 
CF78F7D01 72174381 Be ' > 28.5 ae “XCF 78F 70172174381 3 128*ln2-2%2-49 
DCC858B6 0990429E 3 : ton 5 sG0Ad *XDCC858B60990429E 3 28.5*Ln2 
CF80F7D1 72174031 01 | nase sine “XCF SOF 70172174031 : (high 49 bits of bn2)+2**-49 
FF81898C BOC3ASCA 09 18 ~ “"".QUAD —*XFFB1B9BCBOCSASCA : Un2 = D_.LOG_2_HI 
! DSINHTAB: 
7468 3307 1 WORD *0031727,°0072153 
760F as 1 eWORD %0121466,°0073017 3; DECIMAL: 0.2508223608819151D-07 
EF13 3 1 ORD %*0033470,°016 333 
a H 0 1 «WORD 01544 es ee: 3 DECIMAL: 0.27557298036460860-05 
op 0 3A5 0 1 -WORD *0035120,°0006400 
43A1 D146 1 WORD %0150502,°0041641 : DECIMAL: 0.1984126984813681D-03 
9686 D 4 1 «WORD *003641 01 0ee] 
66A 88 4 1 WORD “01 sel “00 232 3 DECIMAL: 0.8333333333332327D-02 
AAAA 3F 2A 04 1 WORD *0037452,°012525 | 


ete SINH 


88 


00 


Doub 
BECLARA 


B tp 


000006 


f10N 
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f 
5 


10N 


Perr + 4 


une unr— 


ati 


os 
eo 


clarat 


DSINHLEN 


H li 
izperbol ic 


«WORD 
WORD 
- WORD 


—_—_—_—_—_—. 


15 
aSirsyrRain TEEPCHBEE Q4:SD:22 YAMS 82558 46590 sun, 


"Y eogertuteren 


= .> DSINHTAB/8 


; DECIMAL: 0.1666666666666667D+00 
; DECIMAL: 0.D0 
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-SBTTL MTHSDSINH = Standard Double Precision Floating DSINH 


: FUNCTIONAL DESCRIPTION: 

DSINH = double precision floating point function 
DSINH(X) is computed as: 

; ET € geeng 


CALLING SEQUENCE: 
DSINH.wd.v = MTHSDSINH(x.rd.r) 
INPUT PARAMETERS: 


i 

t 

8 

9 

0 

4 

‘ 

8 

9 

0 

4 

5 

6 

7 

8 

1 LONG = 4 ; define get multiplier 
¢ x = 1 * LONG 3; Contents of x is the argument 
4 
5 
6 
7 
8 
9 
0 
4 
5 
6 
7 
8 
9 
0 
4 
5 
6 
7 
8 
9 
0 


00000004 
00000004 


; IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: double precision floating hyperbolic sine of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: 
Signal: MTHS _FLOOVEMAT if 128*ln2 =< iX! with reserved operand in RO/R1 
: (copied to the signal mechanism vector CHF$L_MCH_RO/R1 by LIBSSIGNAL). 
; Associated message is: ‘FLOATING OVERFLOW IN MATA LIBRARY’. Result is reserved 
; operand -0.0 unless a user supplied (or any) error handler changes CHFSL_MCH_RO/Ri 


; NOTE: This procedure disables floating point underflow, enables integer 
overflow. 


= 
oO 
nm 
oO 


ENTRY MTHSDSINH, “M<IV, R2, R3, R4, RS, Rb, R7> 

; standard call-by-reference entry 

; disable DV (and FU), enable IV : 
MTHSFLAG_ JACKET ; flag that this is a jacket procedure in 
6D Q0000000'GF 9 MOVAB G*MTHSSJACKET_HND, (FP) 
pot Seamer address to jacket 
ancier 
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1-008 MTHSDSINH = Standard Double Precision F eet $1355 :88 MTHRTL.SRCJMTHDSINH.MAR; 1 . (4) 
061 
061 ; case of an error in routine 
061 ; If an error, convert signal to user PC 
061 ; and resigna 
56 = 04 BC 70 0061 MOVD @value(AP), R6 3; R6/R7 = {X! = @value(AP) 
0 6 7D 0065 MOVQ R6, RO 
50 8000 8F AA 0068 BICW2 #*x8000, RO 3 RO/R1 = IX! 
3 F 650 «Bl 060 CMPW RO #° x56 80 : compare {X! with 0.25 
1A 18 bre BGEQ  GEG_T0_0.2 > branch if x! >= 0.25 
0074 H 
Baee s iki < 0.25 
74 $ 
0074 
3280 8F 50 81 0074 CMPW RO, #*x3280 3 compare {X! with 2**-28 
06 18 9079 BGEQ GEG_TO_2M28 ; branch if ix! >= 2ee-28 
0078 H 
0078 : IML < 20-28 
0078 3 
0078 
50 56 7D 0078 MOVQ R6, RO ; RO/R1 = X 
04 BO ee RET ; return with result = argument 
007F : 
0077 3 208-28 =< IK! < 0.25 
007F : 
007F 
007F GEQ_TO_2m28: 
50 50 64 OO7F MULD RO,RO ; Get ARG**2 for POLYD 
Al AF 05 50 POLYD RO. #DSINHLEN-1, DSINHTAB 


RO/R1 = SUM(Ci*Xe*i), with 


PVPVPVSVSVIVS B® BB BS BY BE BE AANA AINA PNNININPINNNDY 2 2 OO SS HM DOOOOOOOO 
AALS AN (OOO NAUE WN SO ODNAUE WN 0 ODNOAU EWN CO OONOU EW OO ESENOU SW 
. 


~~ 
w 
oO 
So 
oo 
RIPIPPIAIM INP VUDINININININIPUPYINININIPYINININININPIPININYNIPINININIPYNINININIPIPYIPINYPPYASPIPOPIPUNNPoNPoPy 


0087 3; Last coefficient zero 

50 56 64 0087 MULD R6,RO ; MULD by ARG, and then 

50 56 60 OQO8A ADDD R6,R0 ; Add in ARG with overhang. 

04 008D RET 3 return with result in RO/R1 
OO8E 
O0O8E : 
OO8E 3 6.25 =< ix! 
OO8E : 
OO8E 
OO8E GEQ_T0_0.25: 
FF6D CF 50 71 008 CMPD RO, D_127_L0G_2 3 compare iXi with 127*Ln2 
3B 14 009 BGTR  GTR_TRAN_T27_[0G_2 > branch if ix! > 127#Ln2 
095 3 
0095 3 0.25 =< IX! < 127#ln2 
333 
FF76 CF =650)0=— 71s 0095 CMPD RO, D_28.5 L0G_2 ; Compare iXi with 28.5*\n2, if 
20 14 OO9A BGTR ONLY_ONE_TERM 3 greater, only one of EXP(X) and 

009C ; EXP(-X) is significant. 
009C 3 
009 3 0.25 =< iki < 28.5eLn2 

50 56 7D aes MOVQ R6, RO 3; RO/R1 = X 

OOOO0000"EF 16 O09F JSB MTHSDEXP_R6 3 RO/R1 = DEXP(X) 


ee 
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008 MTHSDSINH = Standard Double Precision F rae yi 94:33:28 EMTHRTL. SREIMTHDS INK MAR: 1 ° (4) 
7E 50 7D OQOAS 57 mMOVQ RO, -(SP) : push DEXP(X) on stack 
5004.8 72 OOK 358 MNEGD avalue(AP), RO i Rove ye 
90000000" © 1g AC 9 JSB MTHSDEXP_R6 + RO/R1 = DEXP(-x) 
2 8 0 6 08 60 SUBD3 RO, (SP)F, R 3: RO/R1 = DEXP(X) = DEXP(-x) 
0 0080 8 A2 0086 61 SUBW «= #*X00B0, RO : RO/R1 = (DEXP(X)=DEXP(-x))/2 
04 0088 66 RET : return with result in RO/R1 
0O0BC 6 
00BC 64 ; 
00BC 65 ; 28.5*ln2 =< IK! <127#Ln2 
Boer 66 ; 
SBE 8% ONLY_ONE_TERM 
00000000" EF 1 pet 63 “ JsB MTHSDEXP_R6 ; RO/R1 = DEXPCIX!) 
a Be ee pla Piette takers 
: negative change sign o 
50 50 72 onc? 372 MNEGD RO. RO : DEXPCTK!) wel Ste 
OOCA 73 POSITIVE: 
50 0080 8F A2 OQOCA 74 SUBW #*x0080, RO 3; RO/R1 = sign(X)*DEXPCIX!)/2 
04 OOCF 75 RET 
0000 276; 
00D0 77 = 127*ln2 < ix! 
0000 78 ; 
0000 79 
00D0 sey GTR_THAN_127_L0G_2: 
FF33 CF 50 71 O0D0 81 CMPD RO, D_128_L0G_2 : Compare {X! to 128*\n2 
1D 18 0005 282 BGEQ ERROR : if 128*ln2 =<ix!, overflow occurs 
0007 283 
0007 284 
0007 285; 
00D7 286 ; 127*Ln2 < ix! < 128*Ln2 
0007 287 ; 
0007 288 
FF3D CF 62 00D7 289 SUBD D_LOG_2. HI, RO ; RO/R1=!Xi-D_LOG 2 HI 
OOOOOO00'EF 16 O0DC 290 JSB MTHSDEXP_R6 ; RO/R1 = DEXPC!XT-B_LOG_2_HI) 
52 50 FF3A CF 65 O0E2 34! MULD3 D_LOG 2_Co, RO, R2 g Re /RS=DEXEC:Xi~D L0G_¢ WT) 9D_.06_2_L0 
50 52 $¢ 00E8 8 SUBD RZ, R 3: RO/R1I=DEXPCiXi-Ln2) 
04 BC 7 OOEB 29 TSTD @value (AP) 3; test the sign of X 
03 18 OOEE 294 BGEQ 10$ 3: branch if X >= 0 
50 50 72 OOFO 295 MNEGD RO, RO ; RO/R1 = si n(x) * DEXPC i Xi =Ln2) 
04 OOFS 296 10%: RET 3 return with result in RO/R1 
OOF4 297 
OOFS 298 ; 
OOF4 99 ; 128*ln2 =< iXi, error 
oor, 809 | 
7E OO°8F 9A OOFS 8 ERROR: MOVZBL #MTHSK_FLOOVEMAT, -(SP) ; condition value 
50 01 OF 79 QOF8 30 ASHQ =#15, #T, RO ; RO/R1 = result = reserved operand -0.0 
OOFC 04 ; ores to signal mechanism vector 
OOFC 05 ; UCHFSL_MCH_RO/R1) so error handler 
OOFC 06 ; can modify the result. 
00000000'GF 01 #4=FB OOFC 07 CALLS #1, G*MTHSSSIGNAL : signal error and use real user's PC 
0103 08 ; independent of CALL vs JS 
04 Ht 1 RET 3 return - RO/R1 restored from CHFSL_MCH_RO/ 
0104 11 
104 \f 
0104 1 END 


MTHSDSINH 
Symbol table 


DSINHLEN 
DSINHTAB 


GEQ_T0_0.25 
GEG"TO-2m08 
GTR-THAN_127_LOG_2 


L 
MTHSSJACKET_HND 
MTHS$S1GNAL 


Initialization 
Command processing 


Pass 
Symbol table sort 
Pass 2 


There were 10 


Page faults 


Symbol table output 
Psect synopsis output 
Cross-reference output 
Assembler run totals 


373 source Lines were read in Pass 
1 page of virtual memory was used to define 1 macro. 


E 16 
; Double Floating Hyperbolic Sine routin Yo i St By 913g) :28 ae 
6-SEP=1984 11:22:46 CMTHRTL 


teeeeeee 
teerenee 
etereree 
00000058 RG 
teetrene 
QOOOOOBC R 
QOOOO0CA R 

= 00000004 


=< 
SOoooooo COOOCOOCOOCOCOCOO 


~~ O —O0— —_—-> 


Allocation 
00000000 < 0.) 
00000104 ( 260.) 


Suceceseusuenesescacesoes$ 


CPU Time 
35 00 :00:00.1¢ 
120 00:00:00.6 
88 00:00:00.90 
0 90:00: 00.01 
68 00:00:00.71 
4 00:00:00.03 
3 BF BO 00' OF 
319 tae 


The working set Limit was 900 pages. A 
4682 bytes (10 pages) of virtual memory were used to buffer the intermediate code. 
pages of symbol table space allocated to hold 22 non-local and 1 local symbols. 

object records in Pass 2. 


+ producing 


00:00:01.7. 
:00:05. 


oo 
oo 
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=~ S5505¢ 
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oo 
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PSECT No. Attributes 
00 ¢ 0.) NOPIC USR CON ABS 
Ot € Bes PIC USR CON” REL 


LCL NOSHR NOE 
LCL = SHR 


Macro V04-00 Page 
« SRCJMTHDSINH.MAR; 1 


XE NORD NOWRT NOVEC BYTE 


EXE RD NOWRT NOVEC LONG 


Brweys lie 
VAX-11 Macro Run Statistics 
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boot eesee rece eee nacee sew em as} 


! ; Macro Library statistics ! 


oe anaunme} 


Macro Library name Macros _defined 
“S25SSDUAZB:(SYSLIBISTARLET.MLB:2 0 itt 0 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHDSINH/OBJ=OBJ$:MTHDSINH MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) #MSRC 
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